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Automated analysis of Fermi-LAT data to trigger ground- based gamma-ray observations 
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Abstract: The Fermi Gamma-ray Space Telescope regularly surveys the entire sky in the energy range between 0.3 and 
100 GeV with an homogeneous coverage. This makes Fermi a very useful guide for ground-based Cherenkov-telescope 
arrays like VERITAS, that are sensitive at energies above ~100 GeV. The VERITAS collaboration uses information from 
Fermi-hAT to select potential targets and identify flaring objects in the GeV band. We describe three different analysis 
pipelines that automatically process and analyze Fermi-hAT data on a daily basis to look for flaring objects: daily light 
curves of selected targets, an all-sky analysis to identify new flaring sources, and an algorithm that looks for clusters of 
high-energy photons in time and arrival direction. 
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1 Introduction 

The Fermi Gamma-ray Space Telescope operates in survey 
mode giving a view of the entire sky in the energy range be- 
tween 100 Me V and 300 GeV every 3h. The sensitivity of 
the Large Area Telescope (LAT) onboard the Fermi satel- 
lite is good enough to detect bright 7-ray sources in time 
scales as short as one day. This is potentially interesting for 
ground-based 7-ray Cherenkov telescopes like VERITAS, 
that can use the all-sky monitoring capabilities of Fermi to 
identify flaring 7-ray sources in the GeV range that could 
potentially be detected at TeV energies from ground. 

The Fen«/-LAT collaboration notifies the community about 
flaring 7-ray sources through Astronomer's Telegrams, 
blogs imiU, and publicly available light curves for a num- 
ber of objects of interest |l3]. However, the information is 
typically made public 1-3 days after the main flaring event 
happened, and sometimes the high gamma-ray activity is 
over. To reduce the latency time between the gamma-ray 
flare and potential VERITAS observations, and also to ex- 
tend the Fermi monitoring to other objects that are of in- 
terest for TeV observatories, the VERITAS collaboration 
has developed different monitoring algorithms that are de- 
scribed in the following section. 

2 Data Analysis 

Three different analysis chains have been developed in- 
side the VERITAS collaboration to monitor ferwi-LAT 
data and promptly identify GeV flaring sources: daily light 
curves of selected blazars, an all-sky analysis to iden- 



tify new flaring sources, and an algorithm that looks for 
clusters of high-energy photons in time and arrival direc- 
tion. The three analysis pipelines automatically download 
Fermi-hPS data and process it using the Fermi Science 
Tools v9rl8 14|. 

2.1 Blazar light curves 

The main goal of this analysis is to produce daily and 
weekly light curves for a list of blazars of interest in or- 
der to identify high flux states in the GeV band. 

The software consists of a series of Python scripts which 
are executed with two separate Perl wrappers. The analysis 
chain consists of two main branches, described below. 

The first branch of the analysis chain performs the time- 
averaged spectral analysis for a given source over an en- 
ergy range between 200 Me V to the highest photon energy 
associated with the source in question. For sources with 
low galactic latitudes, the photon data files (after data se- 
lection) are significantly larger than those of high galactic 
latitude sources. The memory usage for the unbinned like- 
lihood analysis goes linearly with the number of photons. 
As such, it is not possible, for some low-galactic-latitude 
sources, to perform an unbinned likelihood analysis. If a 
photon data file has a size > 50 MB after event selection, 
the time-averaged spectral analysis is performed using a 
binned likelihood analysis rather than the unbinned analy- 
sis. Data within 15° of the source are selected for the time- 
averaged spectral analysis. The source model is built from 
the Fermi IFGL Catalog and includes sources out to an an- 
gular separation of 20°. Spectral parameters for sources 
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beyond 15° are fixed to their IFGL values. All model pa- 
rameters within 15° of the source of interest are initially 
left free. Then the model is fit to the data to find the gen- 
eral location of the minimum in model parameter space. 
All model parameters beyond 10° are then fixed to their 
best fit values. The data is then refit to determine the best 
fit model parameters, and associated errors, for the source 
of interest. The second branch of the analysis chain cal- 
culates the daily and weekly lightcurves for each source. 
Lightcurves of the integral flux above 300 MeV and 1 Ge V 
are produced. The spectral index of the source of interest 
is fixed to the time-averaged value and only the source flux 
is fit. If a particular data point in a lightcurve has a TS 
value less than 0.1, its flux and associated error are set to 

to improve the clarity of the lightcurves. Upper limits 
are not calculated here since they are not used for trigger- 
ing flares. The fitting procedure is the same as that for the 
time-averaged analysis, with one additional step. While 
fitting over smaller time intervals, there will be sources in 
the model that have no associated photons. This leads to 
an excessive number of degrees of freedom that can pre- 
vent the fit from converging. To remedy this, after the fit- 
ting procedure is performed, all sources with TS values less 
than 4 are removed from the model. The data is then refit. 
This process is repeated until there are no longer sources 
with TS values less than the specified threshold. It typi- 
cally does not take more than two iterations to remove all 
non-significant sources. 

To evaluate the significance of a particular flare, the signif- 
icance iS) of a fluctuation in flux is evaluated against the 
average 7-ray flux from the source following the equation: 

S^^l^£^ (1) 

as shown in Figure [T] 

2.2 All-sky monitoring 

This analysis is based on monitoring the all-sky Fermi- 
LAT data daily to watch for increased 7-ray activity above 

1 GeV from any sky direction. This is useful to detect 7- 
ray flares from unknown objects, particularly blazars and 
galactic transients. 

For this analysis, an all-sky exposure corrected map (1° x 
1° bins) is produced including all photons with E > 1 GeV 
collected by Fermi-hP^ during the first 21 months of the 
mission. This map is taken as a static/averaged picture of 
the sky in 7 rays. An equivalent all-sky map is created ev- 
ery day integrating data from the last 3 days. A variability 
parameter (V) is then calculated on a bin by bin basis:. 

^ ^^day 4^21month 

4'21month 

An example of the resulting maps can be seen in Fig |2] 
High values of V are seen as hotspots in the map, and indi- 
cate an increased 7-ray activity in a particular sky direction 



compared to the 21 -month average flux recorded by Fermi. 
In addition to the 3 day integration map, similar maps are 
produced integrating over 1, 7 and 30 days. The differ- 
ent time scales are chosen to test a wide range of variabil- 
ity profiles, going from very fast flares to prolonged high 
states. Whenever a hotspot is identified on a variability 
map, a detailed likelihood analysis of the region is carried 
out, producing the best-fit parameters to the spectrum of 
the flaring source and a light curve to monitor its activity 
above 1 GeV in the last 10 days . This result from the spec- 
tral analysis helps to extrapolate the Fermi derived flux to 
TeV energies, and understand if the source is flaring at lev- 
els detectable by ground-based telescopes. 

2.3 Individual photon counting 

This analysis looks for clusters of photons with E' > 3 GeV 
in the Ferm/-LAT data coming from the direction of a 
known 7-ray source within a short period of time. A cluster 
of high energy photons identifies a source that is in a high 
flux state and naturally selects hard spectrum sources that 
are potentially interesting for VERITAS observations. This 
analysis is particularly suited for sources that show rapid 
variability such as quasars. 

The data processing is quite simple. Photons with E > 
3 GeV are selected. Then the photon list is filtered, 
and only photons within the 95% containment radius 
from the coordinates of the source of interest are kept. 
The containment radius is previously extracted from the 
P6_V3_D1FFUSE instrument response functions |4 | and is 
evaluated separately for front and back converting photons 
as a function of the photon energy (Fig [3]). The remain- 
ing photons, which are coming from the 7-ray source of 
interest with almost negligible background, are plotted as 
a function of arrival time and photon energy (Fig.|4]i. The 
plots show the arrival time and energy (orange dots) of high 
energy photons for the selected sources last 15 days. Alerts 
are typically sent to VERITAS observers when several pho- 
tons are recorded from a particular source within the last 
two days. 

3 Summary 

Three automatic analysis pipelines for Ferm/-LAT data 
have been presented. The different techniques allow the 
identification of 7-ray flares from known sources with high 
sensitivity, the detection of new flaring sources in any sky 
direction, and a prompt reaction to short-lived 7-ray flares 
from rapidly variable objects. The analysis presented here 
allows VERITAS to react in the shortest possible time to 
7-ray flares detected by Fermi, maximizing the chances of 
a simultaneous detection at TeV energies. 
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Figure 1: (top) Gamma-ray light curve for the quasar CGRaBS J1312+4828 in 1-week time bins, (bottom) Significance 
of the fluctuations of each flux point with respect to the source average flux. A dashed (solid) Une indicates 2cr (3a) 
fluctuations. 
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Figure 3: Evaluated Fermi PSF for front and back convert- 
ing photons as a function of energy. 
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Figure 2: {left) Fermi variability map from 2010 December 29th showing the high flux states at i? > 1 GeV from 3C 454.3 
and B2 0200+30, which was observed by VERITAS after approval of a directors discretionary time request, {right) Daily 
light curve for E > 300 MeV measured by Fermi covering the 7-ray flare. 
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Figure 4: Plot showing the arrival time and energy of photons coming from the direction of the quasar 4C 21.35. The 
blue box indicates a tentative trigger criterion. The cluster of high energy photons (orange dots) at MJD 55605 can be 
identified as high 7-ray activity coming from the source. 



